KEK-TH-1080 



O 
O 

(N 
£ 

oo 

(N 
> 

m 

(N 
m 
O 

O 

a: 
43 



Determination of polarized parton distribution functions 
with recent data on polarization asymmetries 

M. Hirai/'Q S. Kumano, 1 - 2 '!] and N. Saito 3 >0 

(Asymmetry Analysis Collaboration) 

'institute of Particle and Nuclear Studies 
High Energy Accelerator Research Organization (KEK) 
1-1, Ooho, Tsukuba, Ibaraki, 305-0801, Japan 
2 Department of Particle and Nuclear Studies 
The Graduate University for Advanced Studies 
1-1, Ooho, Tsukuba, Ibaraki, 305-0801, Japan 
3 Department of Physics, Kyoto University, Kyoto, 606-8502, Japan 
(Dated: March 25, 2006) 

Global analysis has been performed within the next-to-leading order in Quantum Chromodynam- 
ics (QCD) to determine polarized parton distributions with new experimental data in spin asymme- 
tries. The new data set includes JLab, Hermes, and COMPASS measurements on spin asymmetry 
Ai for the neutron and deuteron in lepton scattering. Our new analysis also utilizes the double-spin 
asymmetry for tt° production in polarized pp collisions, A\ L , measured by the Phenix collabora- 
tion. Because of these new data, uncertainties of the polarized PDFs are reduced. In particular, the 
JLab, Hermes, and COMPASS measurements are valuable for determining Ad v (x) at large x and 
Aq(x) at x ~ 0.1. The Phenix ty° data significantly reduce the uncertainty of Ag(x). Furthermore, 
we discuss a possible constraint on Ag(x) at large x by using the Hermes data on gf in comparison 
with the Compass ones at x ~ 0.05. 

PACS numbers: 13.60.Hb,13.88.+e 



I. INTRODUCTION 

The nucleon spin structure has been investigated 
mainly by polarized lepton-nucleon scattering experi- 
ments. After many years of theoretical and experimental 
efforts, our knowledge on the spin structure of the nu- 
cleon is much improved especially in the valence-quark 
sector. However, we still miss a decent measurement of 
orbital angular momenta, transversity distributions, and 
helicity distributions of sea-quarks and gluons in the pro- 
ton. In particular, the polarized gluon distribution has 
large uncertainty. For example, the AAC analysis indi- 
cated that the first moment of Ag(x) is Ag = 0.50 ± 1.27 
at Q 2 ~l GeV 2 0,0 • This represents the largest ambi- 
guity in understanding the longitudinal spin-structure of 
the proton; Ag, the gluon spin contribution to the proton 
can be either positive or negative. 

Fortunately, various experiments started to produce 
data which are sensitive to the gluon polarization. First, 
there were reports on Ag(x) /q(x) from Hermes 3], SMC 
(Spin Muon Collaboration) Q, and Compass collab- 
orations by using high-p^ hadron production measure- 
ments in polarized lepton-nucleon scattering. There are 
also preliminary results [f| from the COMPASS by using 
the data for charmed-meson and high-px-hadron produc- 
tions. Second, the RHIC (Relativistic Heavy Ion Col- 
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lider) Spin experiments started to produce data which 
could constrain the gluon polarization. They arc spin- 
asymmetry measurements in pion 0, and jet Q pro- 
ductions in polarized pp reactions by Phenix and Star 
collaborations, respectively. Because of these new data, 
we expect much improved determination of the polarized 
gluon distribution. 

There are also new measurements on the spin asym- 
metry A\ in lepton-nucleon scattering by JLab (Thomas 
Jefferson National Accelerator Facility) Hall A 10] , Her- 
mes [HI > and COMPASS collaborations in addition to 
the data set used in our previous analysis . These data 
should be useful for a better determination of quark and 
anti-quark distributions, and they may also provide a 
constraint on the gluon polarization through scaling vio- 
lation. In addition, the Phenix ir° data can be included 
in the analysis for a better determination of the polarized 
gluon distribution. 

The major updates from the previous analysis |2j is the 
addition of these new data. Using these reaction data, 
we determine the polarized parton distribution functions 
(PDFs). We discuss analysis results by demonstrating 
the importance of each experiment. In particular, new 
points are the discussions on: 

1. impact of the JLab, Hermes, and Compass data 
on the determination of Aq(x) and Aq(x), 

2. impact of Phenix data on the determination of 
A<7(z), 

3. constraint on Ag(x) from the differences between 
Hermes and Compass data on gf at x ~ 0.05. 



This paper is organized as follows. In Sec. [HJ our 
analysis method is described for determining the polar- 
ized PDFs. Analysis results are explained in Sec. II I II 
The paramctrization results are compared with recent 
spin asymmetry data, the optimum polarized PDFs are 
shown, and they are discussed in comparison with other 
PDFs. Then, details of the polarized gluon distribution 
is discussed in Sec. IIVI The results are summarized in 
Sec. M 



II. 



ANALYSIS METHOD 



Spin-asymmetry data are analyzed to obtain the po- 
larized PDFs. In lepton-nucleon deep inelastic scattering 
(DIS), the spin asymmetry A\ is given by unpolarized 
structure function F2, longitudinal-transverse structure 
function ratio R, and polarized structure function g\ : 



A!(x,Q 2 ) 



9i{x,Q 2 ) 
F 2 {x,Q 2 ) 



2x[l + R(x,Q 2 )}, (1) 



where Q 2 is given by the momentum transfer q as Q 2 = 
—q 2 , and x is the Bjorken scaling variable defined by 
x = Q 2 /(2p ■ q) with the nucleon momentum p. The 
structure function g\ is expressed |l-'l): 

1 n/ r 

9i{x,Q 2 ) = o2Z e n ^C q {x,a s )® [Aq t (x,Q 2 ) 
i=i 

+ Aq z (x,Q 2 )} + AC g {x,a s ) ® Ag(x,Q 2 )X, (2) 



where AC q and AC g are polarized coefficient functions, 
and Aqi, Aq~i, and Ag are polarized quark, antiquark, 
and gluon distribution functions, respectively. The sym- 
bol (g> denotes the convolution integral: 



f(x)®g(x) = 



1 dz 



'<"»(!) 



(3) 



In the similar way, the structure function F% is expressed 
in terms of unpolarized PDFs and coefficient functions. 
The structure functions are calculated in the next-to- 
leading-order (NLO) of the running coupling constant 
a s , and the modified minimal subtraction (MS) scheme 
is used. 

The longitudinal double spin asymmetry for pion pro- 
duction in polarized pp reaction is given by 0, 



.4 



[dcr H 



da^ ] / dpT dAa / dpr 



LL 



[da 



++ 



dcr-| ]/dpx da/dpT 



(4) 



where pr is the transverse momentum of the pion, and 
(Jh 1 h 2 is the cross section with the proton helicities hi 
and h.2- The cross sections Ac and a are defined by 
Act = (a ++ - o+_)/2 and a = {<j++ + o+_)/2. The 



polarized cross sections arc expressed in terms of the po- 
larized PDFs A/ and fragmentation functions [l5| : 



dAa 
dpr 



E 



r/max r 1 

di] / dx a Af a (x a ) 



x / dx b Af b (x b ) 



d(i, z c ) dAa a+b ^ c+x n7r 



d(p T , rj) 



dt 



D?{z c ), (5) 



where d(i, z c )/d(pT, rj) is the Jacobian. The cross section 
dAa a +b^c+x /dt is for the parton-level subprocess a + 
b — ► c + X, and its expression is, for example, found 
in Ref. [T^ . The variables x a and x b are momentum 
fractions for the partons a and b, respectively, and z c is 
the momentum fraction given by z c = p^/Pc with the 
pion momentum p„ and the parton momentum p c . The 
pseudorapidity 77 and t are defined by 



In 



tan — 
2 



t = (Pa - Pc) 2 



(6) 



where 9^ is the angle from the beam direction in the 
laboratory frame. The bounds of the integral over 77, 
Vmax and rjmin , are given by the experimental condition. 
Using Madelstam variables s — (pa +Pb) 2 , t = (pa — 
Pit) 2 , and u = {pb — P-k) 2 , we define variables x\ and X2 



X\ 



u 
s 



XT 



. + '1 



X2 



XT 



(7) 



where xt is given by xt = %Pt / y/s. The lower bounds 
of the integrals and the variable z c are then expressed: 



x a x 2 



1 — x 2 



Xl 



^ + ^. (8) 

Xa X b 



The polarized cross section in Eq. (J5J is calculated 
by these expressions together with polarized PDFs and 
fragmentation functions. The unpolarized cross section 
da/dpT is calculated in the similar way. In calculating 
the NLO cross sections, if-factors are simply multiplied 
in order to obtain them from the leading-order (LO) ones 
as was done in unpolarized PDF analyses [is|. We em- 
ployed 1^=1.0 and 1.6 for the polarized and unpolarized 
cross sections, respectively, basing on more detailed anal- 
yses in Ref. [lq . This approximation inevitably intro- 
duces systematic uncertainties from the pt dependence of 
the K factor and the choice of the PDF set. Such an er- 
ror is estimated to be ^20% in the cross sections, which 
is comparable magnitude to a typical scale uncertainty 
in this energy region. A possible criticism to if-factor 
approach in PDF analysis might be a change in corre- 
sponding x-region due to higher-order corrections, which 
should be the origin of the if-factor. Since the cross sec- 
tion drops rapidly with increasing pt, which is a good 
measure of relevant x region especially for the Phenix 
acceptance covering only central rapidity, changes in the 
cross section would correspond to small change in pt even 
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if the change is fully attributed to a shift in px- For exam- 
ple, 50% changes in the cross section can be obtained by 
10% change in px- Therefore, we think the uncertainty 
involved in the x-determination due to if-factor approach 
is smaller compared to other uncertainties. Clearly more 
precise procedure has to be applied as the precision of 
experimental data improves. 

The polarized PDFs are expressed by a number of pa- 
rameters and unpolarized PDFs at a fixed Q 2 (= Qq): 



by 



Af(x, Ql) = [6x» - k(x v - x»)]f(x, Ql) 



(9) 



where 5, K, v, and \x are parameters to be determined 
by a x 2 analysis, and f(x) is the corresponding unpo- 
larized PDF. Assuming flavor-symmetric polarized anti- 
quark distributions Am = Ad = As = Aq at Qq, we 
determine four polarized distributions, Au v , Ad v , Aq, 
and Ag, which are expressed by the parameters in Eq. 
©■ As usual in this kind of analysis, the first moments 
of Au v and Ad v are fixed by semileptonic decay data 
of octet baryons: J dxAu v — 0.926, J dxAd v — —0.341 

The positivity condition |A/(x)| < f(x) is imposed 
in the initial PDFs. It should be satisfied in the LO, 
whereas it does not have to be strictly satisfied in the 
NLO due to NLO corrections. The details of the NLO 
corrections on the positivity condition are discussed in 
Ref. |l9j ■ We have imposed the condition only for a prac- 
tical reason to avoid unphysical solutions particularly at 
large x. When the precision of experimental data is im- 
proved, we should be able to release this condition. A 
recent unpolarized gluon distribution 0, |2(j becomes 
negative at small x (< 0.01) and Q 2 (~1 GeV 2 ), a spe- 
cial attention should be paid to the condition. In such a 
region, few polarized data are used in our analysis at this 
stage, so that the positivity condition does not affect the 
analysis significantly. 

The parameters of the polarized PDFs are determined 
so as to minimize the total x 2 : 



E 



[^ ata (x,< 



)~A^(x,Q 2 )} 2 



[AA^ ta {x,Q 2 )] 2 



(10) 



where Af ata indicates the spin-asymmetry data for A\ 
and A\ L ; A^ alc is the theoretically calculated asymmetry 
at the same Q 2 point. The polarized and unpolarized 
PDFs at Qg are evolved to the experimental Q 2 point 
by the standard DGLAP (Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi) evolution equations |2l|. The experi- 
mental error AAf ata is given by systematic and statistical 
errors: (AAf ata ) 2 = (AAf at ) 2 + (AAf st ) 2 . The totaly 2 
is minimized by the CERN program library MINUIT [22j . 

The uncertainties of the polarized PDFs are estimated 
by the Hessian method [23j. We express the parameters 
as di (i=l, 2, N), where N is the number of the 
parameters, and the minimum point of \ 2 is denoted a. 
Then, the uncertainty of a distribution F(x) is calculated 



[SF(x)f 



H; 



_! f dF(x,&) 
< v da, 



where is the Hessian. The value of A\ 2 is given by 
Ax 2 = 12.65 so that the uncertainty indicates the one-cr- 
error range for eleven free parameters . The details of 
this A% 2 choice are discussed in a number of publications 
0. 



III. 



RESULTS 



We determine the parameters by minimizing the total 
X 2 in Eq. (|10|) . In order to find an impact of the RHIC- 
A* LL data on the polarized PDF determination, we tried 
two types of analyses: 

• Type 1: with the DIS-^i and RHIC-A£° L data, 

• Type 2: with only the DIS-^4i data. 

For the theoretical ir° production calculations, the kine- 
matical conditions, -^=200 GeV and \r/\ <0.35, are used. 
The initial Q 2 point is taken as Qq = 1 GeV 2 . We 
have tried another fit with Qq = 0.8 GeV 2 , but there 
was no difference in the results. The number of flavor is 
three. The longitudinal-transverse ratio R in Eq. Q is 
taken from the SLAC parametrization [25|. The GRV98 
parametrization j2|j is used for the unpolarized PDFs. 
We also tested other distributions; however, the obtained 
polarized PDFs did not change conspicuously. It is be- 
cause the differences in the PDFs could be absorbed into 
the parametrized function in Eq. by minor adjust- 
ments of the parameters. 

The obtained parameters are tabulated in Tabled ■45| 
Their values are similar to those of the previous version 
. The small- a: behavior of the polarized antiquark dis- 
tributions cannot be determined by the current experi- 
mental data, so that the parameter fiq, which determined 
the small-x behavior, is fixed at [iq = 1 as in Ref. 0. Be- 
cause the first moments of Au v and Ad v are fixed, there 
are no errors in the parameters n Uv and Kd v ■ 

The x 2 values for the used data set are listed in Ta- 
ble [n] The data set includes the measurements from 
EMC (European Muon Collaboration) [27|, SMC pi, 
SLAC-E130, E143, E154, E155 29], Hermes 111 H3, 
JLab Hall- A, [ljj], Compass [l3, and Phenix [ifcollab- 
orations. The additional data, which were not used in 
the AAC03 analysis, are the ones from JLab ^(|> Her- 
mes [H, Compass H2, and Phenix 0. In two dif- 
ferent analyses, similar x 2 /d.o.f. values are obtained. 
As we will see later, two fit results show almost the 
same Au v (x), and slight changes in the central curves 
of Ad v (x), Aq(x), and Ag(x). The major difference is in 
the uncertainty bands for antiquark and gluon distribu- 
tions. 
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TABLE I: Obtained parameters by the \ 2 analyses. 



distribution 


6 


7/ 




M 


Type 1: DIS and tt u data 
Au v 
Ad v 

Aq 

A 5 


0.959 ± 0.099 
-0.773 ± 0.210 
0.780 ± 0.904 
-1.00 ± 4.39 


0.0 (fixed) 
0.0 (fixed) 
1.014 ± 0.180 
2.74 ± 1.30 


0.588 
-0.478 
-75.9 ± 11.2 
252 ± 139 


1.048 ± 0.266 
1.243 ± 0.561 
1.0 (fixed) 
2.70 ± 2.57 


Type 2: DIS data only 
Au v 
Ad v 

Aq 

A 5 


0.958 ± 0.101 
-0.768 ± 0.204 
0.955 ± 0.951 
-1.00 ± 3.28 


0.0 (fixed) 
0.0 (fixed) 
1.014 ± 0.185 
2.30 ± 0.883 


0.585 
-0.474 
-82.7 ± 13.2 
253 ± 153 


1.056 ± 0.272 
1.218 ± 0.552 
1.0 (fixed) 
2.27 ± 1.75 



TABLE II: Numbers of the A\ data and x 2 values are listed for the two types of analyses. The notations p, n, and d indicate 
proton, neutron, and deuteron, respectively. 



data set 


x range 


Q 2 range 
(GeV 2 ) 


No. of data 


2 

X 

Type 1 

ms-Ax+Af L 


Type 2 
DIS-Ai only 


EMC (p) 


0.015-0.466 


3.50-29.5 


10 


4.84 


4.81 


SMC (p) 


0.004-0.484 


1.14-72.10 


59 


55.72 


55.21 


E130 (p) 


0.19-0.64 


5.32-9.91 


8 


4.74 


4.72 


E143 (p) 


0.027-0.749 


1.17-9.52 


81 


61.56 


60.78 


E155 (p) 


0.015-0.750 


1.22-34.72 


24 


32.41 


32.88 


Hermes (p) 


0.033-0.447 


1.22-9.18 


9 


3.32 


3.31 


SMC (d) 


0.0042-0.483 


1.14-71.76 


65 


57.68 


57.16 


E143 {d) 


0.027-0.749 


1.17-9.52 


81 


88.33 


89.13 


E155 (d) 


0.015-0.750 


1.22-34.79 


24 


18.22 


18.26 


Hermes (d) 


0.033-0.446 


1.22-9.16 


9 


12.62 


11.88 


Compass (d) 


0.051-0.474 


1.18-47.5 


12 


9.18 


8.99 


E142 (n) 


0.035-0.466 


1.1-5.5 


8 


2.90 


2.76 


E154 (n) 


0.017-0.564 


1.21-15.0 


11 


3.13 


3.36 


Hermes (n) 


0.033-0.464 


1.22-5.25 


9 


2.14 


2.18 


J-Lab (n) 


0.33-0.60 


2.71-4.83 


3 


2.50 


2.63 


DIS total 






413 


359.29 


358.05 


Phenix (tt°) 






8 


11.18 




total \ A 






421 


370.47 


358.05 


(xVd-o.f.) 








(0.904) 


(0.891) 



A. Comparison with spin asymmetry 
measurements 

We compare the resulted parametrization with newly 
added experimental data of DIS and ir° production in pp 
collisions. Since the agreement with other data is found 
to be maintained, we do not show here. 

First, the parametrization results of the type-1 are 
compared with inclusive DIS data by JLab, Hermes, and 
Compass on the spin asymmetry A\ in Fig. ^ where ac- 
tual differences between the theoretical asymmetries and 
experimental data (^ata _ ^thco-j are s hown. The theo- 
retical asymmetry Af 100 is calculated at the experimen- 
tal Q 2 point of Af ata . The AAC03 uncertainty bands 
are shown by the solid curves and the current AAC06 
type-1 uncertainty ranges are by the shaded areas. The 
uncertainties are calculated at Q 2 —5 GeV 2 . It is obvious 
from this figure that the Hermes and Compass data on 



A\ and Af significantly reduced (about 50%) the uncer- 
tainties of the asymmetries in the region, x ~ 0.2. On 
the other hand, the JLab data play an important role in 
reducing the uncertainty for the neutron in the region, 
x > 0.2. These data are valuable for a better determina- 
tion of the polarized quark and antiquark distributions 
as shown in Sec. IIII Bl 

Second, the analysis results are compared with the 
Phenix data on the pion-production asymmetry A^ L Q 
in Fig. |21 The theoretical curves and their uncertainties 
are shown for the type 1 and 2 analyses. The hard scale 
is taken as Q 2 = p\, in calculating the PDF and their 
uncertainties, and the KKP (Kniel, Kramer, and Potter) 
fragmentation functions |3lj | are used. It is known from 
more detailed analysis that the choice does not change 
the asymmetries |32|. It is clear from Fig. [2 that the 
Phenix data should significantly improve the uncertain- 
ties in the PDFs, mainly the polarized gluon distribution 
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HERMES (p) 



-0.24 
0.2 T 



< -0.2 
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10i 



HERMES (d) 



COMPASS (d) 



-0.2 
(1.2 



JLab (n) 



-0.2^ 

0.001 



0.01 



0.1 



FIG. 1: (Color online) Differences between the data and 
parametrization results. The dotted curves indicate uncer- 
tainty bands of the AAC03 analysis and the shaded areas are 
the bands of the current analysis. 
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FIG. 2: (Color online) Comparison with Phenix 7t°- 
production asymmetry data A\ L . The solid and dashed 
curves indicates type 1 and 2 analysis results, respectively. 
The uncertainty range for the type 1 is shown by the shaded 
band, and the one for the type 2 is by the dotted curves. 



as we will see later. However, the resulted PDFs are still 
consistent with A^ LL = in the entire pt region, and it 
suggests a need of more precise data. Nevertheless, sig- 
nificantly improved constraints on the gluon distribution 
will be discussed basing on this fit results in Sec. IIII Bl 

In order to illustrate the current precision of the struc- 
ture function g\ for the proton, the experimental data 
are compared with the fit results including uncertainty 
bands in Fig. [3J Here, only the type-1 results are 
shown. Experimental measurements are usually listed 



+ 

o 

a. 



0.1 



0.01 



0.7 1 




100 200 



<T[GeV <i ] 



FIG. 3: (Color online) Current type-1 analysis results are 
compared with g\ data. 



by the spin asymmetry so that the g\ data are cal- 
culated by using the unpolarizcd PDFs of the GRV98 
[2^1 and the SLAC parametrization for R with Eq. 0J. 
We notice that the uncertainty bands are still wide at 
large x (>0.5). The smaller x region (x < 0.007) still re- 
mains unmeasured. Although the kinematical coverage 
of the experimental data are being extended, precision 
data are missing especially in the large- X and extremely 
small- a; regions, especially when we compared with the 
HERA data [33]. Wider kinematical coverage is essen- 
tial in extraction of Ag through a scaling violation. Such 
measurement would be feasible at the proposed polarized 
ep-colliders, e-LIC or eRHIC j^. 



B. Polarized parton distribution functions 

The polarized PDFs obtained by the analyses are 
shown in Fig. QJat Q 2 =l GeV 2 , and they are compared 
with the AAC03 distributions. Their uncertainties are 
also shown. The distributions are almost the same in 
both analyses; however, there are much differences be- 
tween the uncertainty bands. Although the Au v uncer- 
tainty bands are the same, the Ad v determination is im- 
proved at x > 0.2 due to the JLab neutron data. The 
polarized antiquark distribution is significantly improved 
at x ~ 0.1 because of Hermes and Compass data on the 
deuteron. 

The most significant improvement is found for the po- 
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FIG. 4: (Color online) Polarized PDFs obtained by the type- 
1 x 2 analysis are shown at Q 2 — l GeV 2 by the solid curves, 
and their uncertainties are shown by the shaded bands. In 
comparison, the AAC03 PDFs are shown by the dashed curves 
and their uncertainties by the dotted curves. 

larized gluon distribution. As clearly shown Fig. the 
uncertainties in gluon polarization are significantly re- 
duced in the type-1 fit, where the All for pion produc- 
tion is included. The central curve also shifted slightly. 
The AAC03 analysis indicated that a negative gluon po- 
larization is possible; however, we clearly see some pref- 
erence for the positively polarized gluon in the type-1 
analysis. We have found a possibility that Ag(x) at large 
x could be constrained by the gf data from the Hermes 
and COMPASS collaborations. The impact of newly ob- 
tained data will be further discussed in Sec. IIVI 

The first moments of Aq(x), Ag(x), and AS (a:) = 

EifeW + are listed at <9 2=1 Gey2 in Table ITTTI 

together with the ones of the AAC03. Because the quark 
and antiquark distributions are almost the same as the 
AAC03, the first moment of Aq and the quark spin con- 
tent AS are similar in the current and AAC03 analyses. 
The only difference is the gluon first moment. First, the 
uncertainty of Ag is significantly reduced in the type-1 
analysis (1.27 — * 0.32). This is because of the Phenix 
data, which constrain the polarized gluon distribution in 
the region x ~ 0.1. Therefore, the Phenix A\ l mea- 
surements play an important role in the determination 
of Ag(x). However, it should be mentioned that an ef- 
fect from the assumed functional form is not included 
in estimating the uncertainty. Furthermore, there is an 
additional scale uncertainty due to the scale error of 
the Phenix data (40%) originated predominantly in the 
beam polarization uncertainty. This uncertainty is esti- 
mated to be ~20% in Ag(x) and its first moement, while 
the quark secter remains intact. 

The polarized PDFs are compared with other recent 
distributions at Q 2 =l GeV 2 in Fig. The solid, dashed, 
dashed-dot, dotted curves indicate AAC06 (type-1 anal- 
ysis), GRSV (Gliick, Reya, Stratmann, and Vogelsang) 
|3^. BB (Blumlein and Bottcher) [3]}, and LSS (Leader, 



TABLE III: The first moments of the obtained polarized 
PDFs at Q 2 — 1 GeV 2 . The current analysis results are 
compared with those of the previous results (AAC03). The 
AE is the quark spin content. 





Aq 


A 3 


AE 


Type 1 


-0.05 ± 0.01 


0.31 ± 0.32 


0.27 ± 0.07 


Type 2 


-0.06 ± 0.02 


0.47 ± 1.08 


0.25 ± 0.10 


AAC03 


-0.06 ± 0.02 


0.50 ± 1.27 


0.21 ± 0.14 



Sidorov, and Stamenov) |38j | distributions, respectively. 
All the distributions agree well except for the polarized 
gluon distribution. However, the uncertainty band of 
Ag{x) is still huge even if the PHENix-pion data are in- 
cluded in the analysis. The different distributions agree 
each other because they are roughly within the uncertain- 
ties. The various parametrizations have different analy- 
sis methods such as in the used data set, x-dependent 
functional form, treatment of higher-twist effects, and 
positivity condition; however, the obtained distributions 
are similar. There are other recent parametrization stud- 
ies |39| with semi-inclusive data by de Florian, Navarro, 
and Sassot and also in the statistical-parton-distribution 
model by Bourrely, Soffer, and Buccella although they 
are not shown in the figure. There are also other 
parametrization studies in Ref . |40| . 

Next, the polarized antiquark distribution of the type-1 
is compared with experimental antiquark-distribution ra- 
tios Au/u, Ad/d, and As/s by the Hermes collaboration 
[41| in Fig. The Hermes measured semi-inclusive spin 
asymmetries for pion and kaon productions in positron 
and electron DIS with proton and deuteron targets. They 
analyzed the semi-inclusive data together with the in- 
clusive ones on A\ for obtaining the quark and anti- 
quark distribution ratios, Aq(x)/q(x) and Aq(x)/q(x), 
in the LO. Effects of NLO radiative corrections are small 
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FIG. 5: (Color online) Comparison with other polarized 
PDFs at Q 2 = l GeV 2 . The type 1 distributions and their 
uncertainties are shown by the solid curves and bands. The 
others are the GRSV, BB, and LSS parametrizations. 



G 



IV. POLARIZED GLUON DISTRIBUTION 
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FIG. 6: (Color online) Comparison with the Hermes semi- 
inclusive DIS data on the antiquark distribution ratios Au/u, 
Ad/d, and As/s. The solid curves are the type-1 analysis 
results at Q 2 =l GeV 2 and the uncertainty ranges are shown 
by the shaded bands. 



in comparison with the experimental errors because the 
distribution ratios are taken. Therefore, although the 
experimental ratios are obtained by LO analyses and 
the theoretical ratios are calculated in the NLO, they 
can be compared with each other. In the AAC analy- 
ses, flavor symmetric antiquark distributions are assumed 
Au(x) = Ad(x) = As(x) at Q 2 =l GeV 2 ; however, one 
should note that the antiquark distributions are different 
(Au(x) ^ Ad(x) ^ As(x) at Q 2 ^ GeV 2 ) because of Q 2 
evolution effects in the NLO |42j. In the figure, the type-1 
distributions are shown aXQ 2 —l GeV 2 ; however, Q 2 vari- 
ations are small because the ratios are taken. The theo- 
retical antiquark-distribution ratios agree with the Her- 
mes data if the experimental errors and the parametriza- 
tion uncertainties are considered. At this stage, the semi- 
inclusive data are not included in the global analysis. 
The data are fully consistent with our parametrization. 
It is obvious that the data would impose significant con- 
straints when the data precision is improved. 

The AAC00 and AAC03 codes are available at the web 
site 2] for calculating the polarized PDFs. However, such 
a library is not prepared for the current analyses because 
the Phenix data utilized here are still preliminary re- 
sults. In particular, the data will become more accurate 
in the final paper, so that the Ag(x) uncertainty should 
become smaller. After the data are finalized, we plan to 
provide a code for the polarized PDFs. 



In this section, the details are discussed on the deter- 
mination of the polarized gluon distribution. First, wc 
explain the possibility that the recent Hermes and COM- 
PASS measurements on Af could impose a constraint for 
Ag(x) at large x because of their Q 2 differences. Second, 
the possibility of negative gluon distribution (Ag(x) < 0) 
at small x is studied by including the Phenix it data 
in the analysis. Third, the polarized gluon distributions 
obtained by the current analyses are compared with Her- 
mes, SMC, and Compass measurements on the ratio 
Ag(x)/g(x). 



A. Polarized gluon distribution at large x 

We discuss the possibility that precise measurements 
of g\ at x ~ 0.05 could impose a constraint on Ag(x) 
at larger x. As shown in Eq. J2J), the structure func- 
tion is given by the convolution integral of the po- 
larized PDFs with corresponding coefficient functions. 
In particular, we discuss an effect of the gluon term 
(x/ z)Ag(x/ z)ACg(z), where z is the integration variable 
in Eq. ©. 

The parametrization results of the type 1 are compared 
with inclusive DIS data particularly by JLab, Hermes, 
and COMPASS on the spin asymmetry A\. Then, the 
parametrization is consistent with the SLAC-E154 neu- 
tron data and COMPASS deuteron data. Although the 
agreement is excellent also with the Hermes proton data 
for A\, the curve deviates at small x (0.03 < x < 0.07) 
in the Hermes deuteron data (Af). 

It is noteworthy that the central curve of our fit re- 
sults deviates from Hermes Af at small x, although 
a reasonable agreement is obtained with Compass Af 
in the same x-region. The deviation from the Hermes 
data points becomes larger when we artificially remove 




FIG. 7: (Color online) The Hermes and Compass data on 
the spin asymmetry of the deuteron (Af) at small x axe shown 
in comparison with the parametrization results. The solid 
curves indicate the type-1 asymmetries and the dashed ones 
are obtained simply by removing the NLO gluon term ( AC g = 
0) in the structure function g± . 
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FIG. 8: (Color online) The polarized gluon distributions and 
coefficient function are shown at £=0.001, 0.05, and 0.3 for 
the type-1 parametrization. 



the gluon contribution to gf , which is shown as dashed 
curves in Fig. Because the Q 2 values are quite dif- 



ferent: Q 



2 

Hermes 



1 GeV 2 and Q 



Compass 



6 GeV 2 , we 



could attribute the observed discrepancy to the Q 2 dif- 
ference. There are two possible effects of such kind; one 
is the higher-twist effect as pointed out recently in Ref. 
|38|. The other, which we found in this analysis, is due to 
a positive gluon polarization at large x, which we discuss 
below. 

We recall the gluon contribution to the gi is 
(xj z)Ag(x/ z)ACg(z) integrated with respect to z over 
(a;, 1). In this picture, the gluon with the momentum 
fraction x/z would split into the quark with a momentum 
fraction x and anything else to contribute to the gi(x). 
Figure [S] displays the type-1 polarized gluon distribution 
(x/z)Ag(x/z) for selected x- values and the polarized co- 
efficient function AC g {z) as functions of z. The AC g (z) 
is positive in the region, 0.05 < z < 0.7, whereas it is 
negative in other z region. The function (x/z)Ag(x/z) 
for x = 0.001 peaks at around z = 0.002 and it moves to 
a larger z as x becomes larger (x=0.05 and 0.3). There- 
fore, the gluon contribution to the g\ is negative in the 
small- a; region, positive at medium x (~ 0.05), and back 
to negative at large x. 

The Hermes data shows a discrepancy from our 
parametrization result with the gluon contribution elimi- 
nated especially in the region, 0.033< x <0.065. This 
is the x region where the gluonic contribution from 
Ag(x) > is positive. If the Ag is negative, the gluonic 
contribution to will be negative so that the discrep- 
ancy becomes larger. Therefore, we could interpret the 
deviation of Hermes data as an indication of positive 
Ag(x) at large x. 

Obviously further precision studies are necessary to 
identify the source of the discrepancy. Determination of 
the sign of the gluon polarization will be helpful. More 
precise data on the Q 2 dependence of g\ would be cer- 
tainly helpful, too. 



B. Polarized gluon distribution at small x 

We have discussed the polarized gluon distribution at 
large x basing on Q 2 dependence of Af . On the other 
hand, the gluon polarization in small- x region is not ob- 
vious. Most of the global analyses showed Ag(x) > 
in the small- x region, and estimated uncertainty turned 
out to be small including our previous analysis, BB, and 
LSS.@We are going to show that different small- a; behav- 
ior is still allowed within the current framework, and that 
the uncertainty band changes drastically. 

We start from the possibility of negative gluon po- 
larization inspired from the two possible solutions of 
Ag{x)/g{x) to reproduce A\ L measured by Phenix. 
Since the pion production is dominated by gg scattering 
at small pt, the asymmetry A\ L approximately depends 
on quadratic function of Ag(x)/g(x). Indeed, two solu- 
tions, positive and negative Ag(x) / g(x) can be obtained, 
when only A\ L data are used to constrain Ag(x)/g(x) 
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FIG. 9: (Color online) Theoretical asymmetries are shown 
for the type 1, 2, and 3 parametrizations in comparison with 
the Phenix data of A^ L . Only the small-pT region is shown 
in order to see the differences between the theoretical asym- 
metries at small pr- 
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FIG. 10: (Color online) The polarized gluon distributions are 
shown at Q 2 =l GeV 2 . The solid, dashed, and dotted curves 
indicate the type 1, 2, and 3 distributions, respectively. The 
uncertainty ranges are shown by the shaded band, the dotted 
curves, and the thin solid curves for the type 1, 2 and 3, 
respectively. 
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[TH |43| . The Phenix A^' l is sensitive to the gluon 
around x ~ 0.1, and the negative Ag(x)/g(x) solu- 
tion can coexists with the previously mentioned positive 
Ag{x) at large x. 

To explore this possibility, we have performed another 
global fit by assigning the initial distribution, Ag(x) = 
—Ag(x)typ e -i, at Qq without modifying the parametriza- 
tion in Eq. (JHJ). The fit was successful, and there was no 
practical difference in y 2 . The resulted parameter set is 
referred to as type 3. The fit results are shown in Fig. 
El in comparison with the type-1 and 2 parametrizations 
and the Phenix experimental data on the double spin 
asymmetry A\ L . For explaining the negative asymmetry 
at pt=2.38 GeV, a significant change of the functional 
form of Ag{x) is needed in the region, 0.06 < x < 0.2. 
There is no noticeable change in the quark distributions 
and only xAg(x) is shown in Fig. 1101 As can be expected 
from the discussion above, the Ag(x) is negative up to 
x ~ 0.1 and retains positive value towards larger x region 
especially to reproduce Q 2 dependence of Af. 

Three types of the gluon polarization Ag(x) are dis- 
played with their uncertainties in Fig. 1101 The uncer- 
tainty bands are reduced in the analyses with the pion 
data. In principle, Ag(x) should be constrained only 
in the PHENIX kinematical range x ~ 0.1. However, 
the reduction is also found in the region x > 0.2. This 
is because a smooth functional form is assumed in the 
analysis, where the constraint in the x ~ 0.1 region 
also affects the uncertainty estimation in other x regions. 
Because the gluon-gluon subprocesses dominate, the po- 
larized cross section for pion-production depends on Ag 
quadratically. It leads to two independent solutions, the 
type 1 and type 3. In the region x > 0.1 partially cov- 
ered by the PHENIX data, similar uncertainties are ob- 
tained and they cover both solutions. However, the type- 
3 Ag(x) resides outside of the uncertainty bands of type 1 
and type 2 at small x, which reveals that the uncertainty 
depends on the functional behavior of Ag{x). The errors 
due to this functional form are not included in showing 
the uncertainty bands in our analysis and also other ones 
[37l l3^ | , because there is no systematic way to include 
them. 

The large uncertainty in type-3 Ag(x) in the smaller-x 
region is also reflected to the first moment: —0.56 ±2.16. 
The large error comes from the small- a; region, where 
there is no experimental data which constrains the small- 
x behavior of Ag(x). In fact, if the type-3 Ag(x) is 
integrated over the region, 0.1 < x < 1, covered by 
the DIS and pion-production data, the first moment be- 
comes 0.32±0.42 which is comparable to the type-1 value, 
0.30 ± 0.30, in the same x region. Only the Q 2 evolution 
and smoothness of functional forms are the possibilities 
to provide some constraints in the small- a; region. Please 
note that there is no direct experimental constraint for 
Ag(x) at small x, therefore it is highly unconstrained. In 
order to clarify the situation, we need a definitive mea- 
surement of the sign of the gluon polarization around 
x < 0.1. The higher energy run (y^s = 500 GeV) at 



RHIC and eRHIC/e-LIC would play important roles in 
such measurements. 



C. Comparison with high-pT hadron-production 
data in DIS 

There are reports on the polarized gluon distribution 
from high-pT hadron production in DIS by Hermes 0, 
SMC 4], and COMPASS collaborations. Their data 
are shown by the ratio Ag(x)/g(x). There are also pre- 
liminary data by the COMPASS [6J for Ag(x)/g(x) from 
analyses of high-Q 2 (>1 GeV 2 ) hadron-production and 
charmed-hadron-production data. 

We compare our results with those experimental data 
in Fig. The parametrization results for the type 1, 2, 
and 3 are shown together with their uncertainty bands. 
The high-pr hadron data in DIS impose a constraint for 
the polarized gluon distribution in the x range, 0.07 < 
x < 0.2. We find that three fit results are consistent with 
these data within both uncertainties. 

It is interesting to note that the Ag(x) / g(x) distribu- 
tions obtained by the analyses are similar to a theoretical 
prediction by Brodsky and Schmidt |44| in the region, 
z~0.1. 
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FIG. 11: (Color online) The theoretical ratios Ag(x) / g(x) 
of the type 1, 2, and 3 are compared with high-pT hadron- 
production data measured in DIS by Hermes, SMC, and 
COMPASS collaborations. The solid, dashed, and dashed-dot 
curves indicate type 1, 2, and 3 parametrization results, re- 
spectively. The uncertainty ranges are shown by the shaded 
band, the dotted curves, and the thin solid curves for the type 
1, 2 and 3, respectively. 



V. SUMMARY 

We have performed the global analysis of spin asymme- 
try data from DIS and neutral-pion production in polar- 
ized pp scattering. First, we showed two types of analysis 
results by using DIS and Phenix pion data and by us- 
ing only the DIS data. In comparison with the previous 
version (AAC03), new data are added from the JLab A™, 
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Hermes and Compass Af, and Phenix measure- 
ments. We found that the JLab neutron data contribute 
to reduce the uncertainty band of the valence distribu- 
tion Ad Vl the Hermes and Compass data to improve 
the antiquark distribution Aq in the region, x ~ 0.1. 

There are two important consequences on the polar- 
ized gluon distribution. First, the differences between the 
Hermes and Compass Af data could indicate a positive 
gluon polarization at large x because of their Q 2 differ- 
ences. Second, the Phenix A^ ll measurements improve 
the determination of Ag(x) significantly. In fact, the un- 
certainty band of Ag(x) becomes smaller by 60%. 
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